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Quantum gas-liquid condensation in an attractive Bose gas 
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Gas-liquid condensation (GLC) in an attractive Bose gas is studied on the basis of statistical 
mechanics. Using some results in combinatorial mathematics, the following are derived: (1) With 
decreasing temperature, the Bose-statistical coherence grows in the many-body wave function, which 
gives rise to the divergence of the grand partition function prior to Bose-Einstein condensation. It 
is a quantum-mechanical analogue to the GLC in a classical gas (quantum GLC). (2) This GLC is 
triggered by the bosons with zero momentum. Compared with the classical GLC, an incomparably 
weaker attractive force creates it. For the system showing the quantum GLC, we discuss a cold 
helium 4 gas at sufficiently low pressure. 

PACS numbers: 67.20.+k, 67.40.-w, 64.10.+h 



I. INTRODUCTION 

Gas-liquid condensation (GLC) is one of the universal 
phenomena occurring both in the classical and the quan- 
tum world. GLC is regarded as a singularity appearing 
in the pressure- volume curve determined by the equation 
of state: 
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The question is whether the grand partition function 
Zy (/i) shows the singularity [H . In ordinary GLC such 
as vapor condensation to water, the GLC occurs at rel- 
atively high temperature compared with the quantized 
energy. Hence, this GLC belongs to the classical phe- 
nomenon. 

When GLC occurs at low temperature as in helium gas, 
however, the situation is different. One must seriously 
take into account the influence of quantum statistics on 
the GLC. In this paper, we call such a GLC quantum gas- 
liquid condensation to emphasize the key role quantum 
mechanics plays in the phenomenon, and we study its 
peculiar nature. 

A prototype of the quantum GLC is the GLC in an at- 
tractive Bose gas at sufficiently low temperature. For the 
relationship between the GLC and Bose statistics, it has 
been often discussed in connection with Bose-Einstein 
condensation (BEC) [3]. The following points, however, 
should be noted. (1) No BEC gas with an attractive 
interaction remains stable in thermal equilibrium. (The 
velocity of sound propagating through the condensate be- 
comes imaginary (2) Since one of the two phases 
does not exist in the BEC phase, a transition between 
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FIG. 1: The phase diagram of helium 4 at low pressure. (In 
the inset, the wavy lines and the dotted lines around particles 
schematically represent the wave function and the zero-point 
motion, respectively.) 



them is impossible. (3) Rather, the quantum GLC oc- 
curs in the normal phase close to the BEC transition 
point . In this region, the coherent many-body wave func- 
tion is composed of many Bose particles, but it does not 
yet reach a macroscopic size. With decreasing temper- 
ature, the bosons' kinetic energy approaches zero, but 
they do not experience a repulsive force until their dis- 
tance reaches the hard-core radius. When a weak attrac- 
tive force acts on such bosons, its influence is drastic. A 
dilute Bose gas will undergo GLC prior to BEC [5j . Its re- 
markable difference from classical GLC is that, under the 
influence of Bose statistics, a negligibly weak interaction, 
which plays almost no role in the classical phenomenon, 
causes drastic changes of the system. 

For the system showing the quantum GLC, we know 
an old example. At pressure lower than 0.05 atm (triple 
point), a cold helium 4 gas in the normal phase under- 
goes a phase transition directly to a liquid in the super- 
fluid phase as depicted in Fig.l @. At sufficiently low 
pressure, a gas is so rarefied that atoms experience a 
weak attractive interaction with distant atoms. Because 
of this weakness, the system remains in the gas state to 
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such a temperature that its kinetic motion is affected by 
quantum statistics. Hence, it is probable that quantum 
statistics directly determines the character of this GLC 
0. 

An essential point to formulate the quantum GLC is 
that the coherent many-body wave function, which is 
symmetrical with respect to interchanges of a large num- 
ber of particles, grows in Zy(p) at low temperature. At 
the vicinity of Tbec m the normal phase, it will give 
rise to the instability of the system. A precise formu- 
lation of this quantum effect is worth studying closely. 
The peculiar nature of Bose statistics manifests itself in 
momentum space. Hence, if one formulates the pertur- 
bation expansion of Zy{p) in momentum space, one may 
be able to obtain a straightforward expression of the in- 
stability. Along this line of thought, statistical mechanics 
of the quantum GLC was recently considered [§]. 

In this paper, using a different method from that in 
Ref.Q (see Appendix. A), we give a statistical-mechanical 
model of the quantum GLC . Mayer and others con- 
sidered the divergence of the perturbation expansion of 
Zv{p) [H- We will consider the quantum GLC along 
Mayer's viewpoint Using some results in combinato- 
rial mathematics, we will derive a set of upper and lower 
bounds for Zy{\i). By studying the radius of convergence 
for both bounds, we find that the divergence of Zy{p) oc- 
curs prior to BEC in cooling, which leads to the onset of 
the quantum GLC . 

An outline of this paper is as follows. Sec. 2 describes 
the Bose-statistical coherence appearing in Zy (p) at low 
temperature. Sec. 3 develops a perturbation formalism of 
an attractive Bose gas along the viewpoint of Sec. 2, and 
leads to the instability. Sec. 4 compares the classical GLC 
with the quantum GLC , and discusses some differences 
between an attractive Bose and an attractive Fermi gas. 



II. BOSE-STATISTICAL COHERENCE 

A. A model 

We consider a spinless attractive Bose gas with a 
weakly attractive interaction Hn 



H = e p ala p + U ^ a l- 
p p,p',q 
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FIG. 2: A pattern of elementary polygons in coordinate 
space, representing a set of permutations P in the exponent of 
Eq.(8). The size of each polygon reflects a size of the coherent 
many-body wave function in Zo((J,). 



perturbation theory with respect to the interaction Ha 
' (-l) n 



Z v (n) = Z (fi) Y 



ri 



n=0 

x dpi--- I d0 n (TH it (/3 1 )---H u (/3 n )), (4) 
Jo Jo 

where Zq(/j,) is the grand partition function of an ideal 
Bose gas. ((—1)" comes from the Boltzmann factor.) 
The emergence of the Bose-statistical coherence appears 
not only in Z (jj) but also in (TH it (Px) ■ ■ ■ H it (j3 n )}. 



B. The unperturbed part Zo(/j,) 

Figure. 2 schematically shows the coherent many-body 
wave functions in coordinate space. (Thick arrows repre- 
sent permutation operations of Bose particles.) Matsub- 
ara [l(| and Feynman [ll| associated the grand partition 



function Zq([i) = (1 — e 



) 1 with such coher- 



ent wave functions as follows. The partition function 
Zq (N) of N Bose particles is given by the density matrix 
p(xi, . . . x N ;x[, . . . x' N ; (3) as [12 1 



Z (N) 



per 



.Xjy] Pxi, . . . Pxjv; 0)d x 



where P denotes permutations within N particles, and 
p is summed over all possible permutations. For the un- 
perturbed part, the density matrix satisfying the "ideal 
gas" equation 



The situation in a low-density gas, in which particles ex- 
perience a negligibly weak interaction since they are far 
apart, and experience it only when a particle encounters 
other particles, enables us to use a contact interaction 
U . To characterize the initial stage of the instability 
occurring in a rarefied gas, we ignore the short-range re- 
pulsive force which plays an important role after the sys- 
tem reaches the high-density state. We regard ^ p e p aj,a p 
as the unperturbed Hamiltonian, and obtain the grand 
partition function Zy(^)=Tr exp[—(3(H — pN)\ by the 
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One therefore obtains 
Z (N)= ' 



N\ \2n f3h 2 



J E ex P 

per 
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d N Xi . (8) 



This expression of Zq{N) is interpreted as follows. For 
N = 3 and P = (1 -> 2,2 -> 3,3 -> 1), the exponent 
of Eq.(8) includes (x\ — a; 2 ) 2 + (a; 2 — x 3 ) 2 + (x 3 — xi) 2 . 
These permutations represent a triangle in Fig. 2, which 
corresponds to a small coherent wave function. In gen- 
eral, the same particle appears only twice in T,(xi — Pxi) 2 
, first at an initial position Xi and second at one of the 
final positions Pxi. The intermediate permutations be- 
tween the initial and final x s form a closed graph, which 
is visualized as a polygon in Fig. 2. As a result, one 
finds that many polygons are spreading out in coordi- 
nate space. (We call this an elementary polygon.) The 
size of the polygon (the number of its sides) corresponds 
to the number of particles, hence the size of the coher- 
ent wave function obeying Bose statistics {coherence size 
[lH). Each particle belongs to one of these polygons, and 
polygons do not share the same particle. A set of per- 
mutations in Eq.(8) corresponds to a pattern of polygons 
like FigJ?. 

Ref. 1(| 1I| calculated Z (N) by the following geomet- 
rical consideration. A polygon of size s appears in the 
integrand in Eq.(8) as 
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Assume a polygon distribution as {£i, £2, ■ ■ ■ , £,s, ■ ■ ■ }, in 
which elementary polygons of size s appear £ s times, be- 
ing subject to N — £ s s£ s . Consider the number of all 
configurations associated with {£ s }, and denote it with 
B(£i, ...,£ s ,.. .). One can rewrite Eq.(8) as 



Z (N) 



1 

{£»} 
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(1) -B(£i, . . - ,£«>• • •) is estimated as follows [Tl|. As- 
sume N particles in an array. The number of ways of 
partitioning them into {£ s } is given by Nl/H s !; s l. An 
array of s particles corresponds to a polygon of size s. 
For the coherent wave function, it does not matter which 
particle is an initial one in the array (circular permuta- 
tion). Hence, A!/n s ^ s ! must be multiplied by a factor 
1/s for each £ s , with a result that 



Nl 



(ii) 



With Eq.(ll) in Eq.(lO), one obtains 
Z (N) 



e n 



(12) 



where A = {2irf3h 2 /m) 1 / 2 is the thermal wavelength. 

(2) Using Eq.(12), we obtain the grand partition func- 
tion Zo(n) = T,;\[Zo(N)e^ tJ ' N , in which summation over 
N changes to free summation over £ s from to 00. Sub- 
stituting N = Y, e s$,s into \ 3N in Eq.(12) and e 13 ^ yields 



(13) 



The f s in Eq.(9) is estimated with the convolution theo- 
rem (see Appendix. B), with which one obtains 



Z ((i) = exp 



E 



A. 
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(14) 



The first and the second term in the exponent come 
from p = and p bosons, respectively, both of which 
are expansions with respect to the coherence size s [14| . 
Using 



E-=mi 



(15) 



s=l 



and 



E^--^/ y 2 Ml-*e-»)dy, (16) 



Eq.(14) agrees with Z (n) = - el 1 



Using Eq.(14) in Eq.(l), the coherent many-body wave 
function appears in the equation of states as follows: 



P 



lim 



1 ^ h(s) 



(17) 



where h(s) is the size distribution of the coherent wave 
function (1/s on the right-hand side is a geometrical fac- 
tor.). In view of Eq.(14), one has h(s) — e^^ s for p = 
(1°) and h(s) = (A s U/A 3 )(e /3 ' JS /s 1 ' 5 ) for P^- The chemi- 
cal potential p, determines the size distribution. At high 
temperature (/x -C 0), h(s) falls off exponentially with a 
size s, and large polygons are therefore negligible (Boltz- 
mann statistics). With decreasing temperature (/z — > 0), 
the exponential dependence on s weakens as e^ MS — > 1, 
thus making large polygons significant. (At the BEC 
transition point, its size dependence disappears forp = 0, 
and the macroscopic wave function therefore contributes 
to Zq(/j,) as much as the microscopic one.) It should be 
noted that, even in the normal phase, the large but not 
macroscopic coherent many-body wave functions exist at 
the vicinity of the BEC transition point, and make con- 
siderable contributions to the grand partition function. 
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FIG. 3: Some ring diagrams in momentum space, composed 
of bubbles of bosons (an example of n—7 terms in Eq.(4)). 
Dotted lines represent the attractive interaction in Eq.(3). 



FIG. 4: A polygon cluster including a square and a hexagon, 
which are made by exchanging particles between two bubbles 
and between three bubbles in Fig. 3, respectively. It belongs to 
the (2, 1, 1, 0, 0, 0, 0) pattern, an example of K" 1 ■ ■ ■ K^'-type 
graphs. 



C. (THu(/3i) ■ ■ ■ Hit(/3 n )) 



In an attractive Bose gas, we must extend the 
concept of Bose-statistical coherence from Z (^) to 
(TH it ((3i) ■ ■ ■ H it (f3 n )) in Eq.(4). 

At high temperature, relatively simple diagrams are 
significant in (TH it ({3i) ■ ■ ■ H it (f3 n )}. A typical example 
in momentum space is the sum of ring diagrams com- 
posed of bubbles like Fig. 3. A thin solid line with an 
arrow represents a boson propagator, and dotted lines 
denote the attractive interaction. 

Assume that (TH it (l3i) ■ ■ ■ H it (f3 n )} is decomposed into 
hi + h 2 + ■ ■ ■ disconnected diagrams containing closed 
loops, in which each of h\ diagrams contains m\ interac- 
tion lines, each of h 2 diagrams contains interaction 
lines, and so on (n = J2i^i m i)- The total number of 



such diagrams is n\/h\\h2 



A disconnected diagram 



with m interaction lines contributes to Zv{n) as 
(~l) m 



rn 



x f d^--- d(3 m (TH it ((3i) ■ ■ ■ H it ((3 m )). (18) 
Jo Jo 

Hence, a sum of all disconnected diagrams 

hi h 2 

replaces the multiple integrals with (—1)" in Eq.(4), and 
free sums over hi replace the sum over n as 



Z v {p) = Z (fi) exp (Si + ^ 2 + 



(20) 



With decreasing temperature, the wave function in 
(TH it ((3i) ■ ■ ■ H it {j3n)) becomes symmetrical with respect 
to interchanges of particles. The sum of rings in Fig. 3 



changes as follows: When one of two particles in one 
bubble and that of another bubble have the same mo- 
mentum {p — p'), an exchange of these particles by the 
thick arrows yields a square as in Fig. 4, which links two 
small rings. The inclusion of such diagrams in Zv(f^) 
ensures the principle of Bose statistics that many identi- 
cal particles have the tendency to occupy the same state 
(P = P')- 

Similarly, when three particles arc exchanged between 
three bubbles in the lower part of Fig. 3, a hexagon is 
yielded as in Fig. 4. Further, when two particles are in- 
terchanged between the square and the hexagon in the 
cluster, a decagon is yielded. Such a sequence will con- 
tinue to a single large polygon in Fig. 5. Starting from 
n bubbles, n times of particle interchange between two 
bubbles leads to a single large polygon consisting of 2n 
bosons (2n gon) . The size of the polygons (the number of 
its sides) reflects the number of particles, hence the size 
of the coherent wave function. With decreasing temper- 
ature, the large coherent wave-function becomes impor- 
tant not only in Z (/j,) but also in (TH it ([3i) ■ ■ ■ H it ((3 n )}. 

In the unperturbed part Z (fi), we consider the ele- 
mentary polygons. Here, we find another type of poly- 
gon at each order of the perturbation. From now on we 
call this type of polygons an interaction polygon. Both 
polygons are different, complementary expressions of the 
coherent many-body wave function. In contrast with the 
elementary polygons, this type of polygon has an interac- 
tion line at each vertex, thus being connected with each 
other. As depicted in Fig. 4, a 2s-sized polygon is com- 
posed of s bosons with (j>i, U), and another s bosons with 
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which is the analogue of Eq.(lO) for Z (N). H lt in Eq.(3) 
is a function of the distance between particles in coordi- 
nate space. When the integral in coordinate space is 
performed in (THi t (Pi) ■ ■ ■ Hi t ((3 n )), relative coordinates 
of particles are used, thus reducing multiplicity of the 
integral by one, with a further factor V entering in the 
right-hand side of Eq.(22), which represents the transla- 
tion of the polygon cluster . 

For evaluating S n , one must consider the influence of 
Bose statistics on [] s (if s ) l ' > and make a graphical con- 
sideration of D n (i>i, . . . , v s , . . .). 



Effects of Bose statistics 



FIG. 5: A single largest polygon (14-gon) as a final result 
of particle interchanges. It belongs to the (0, 0, 0, 0, 0, 0, 1) 
pattern, an example of A^-type graphs. 



[~{Pi + Si)' — (h + m i)]j with an expression given by 



K„ = 



s i 

/ 



(21) 
\ 



fe^^+'y (e-(i*+«B) - A*) 



. Ti(k + m i) 
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where l,m is an even integer including zero. While the 
elementary polygons in Zq((j,) represent a series of succes- 
sive one-particle permutations like Fig. 2, the interaction 
polygons represent intertwining of many permutations. 



III. FORMALISM 

To formulate the concept of polygons in 
(THit(Pi) ■■■ Ha (Pn)), h is useful to make reference 
to the description below Eq.(9). In a polygon cluster, 
let us assume a distribution of interaction polygons as 
{v s } — {^1,1%...}, in which a polygon K s of size 2s 
appears u s times. (For example, {u s } — {2,1,1,0,...} 
in Fig. 4.) They are connected to each other by n 
interaction lines, in which 2s interaction lines emerge 
from 2s-sized polygon. Hence, n — sv s must be 
satisfied. 

Consider the number of all configurations correspond- 
ing to {v s }, and denote it by D n (y\, . . . ,v a , . . .). One can 
rewrite S„ in Eq.(20) with Eq.(18) in terms of interaction 
polygons as follows: 



S„ = V— V* D n (ui, 
{".} 



)k?...k:° 



(22) 



In Y\ s {K s ) Vs of Eq.(22), Bose statistics manifests itself 
in the following manner. 

(1) With decreasing temperature, the feature of Bose 
statistics, accumulation of particles in the low-energy 
state, becomes apparent, thus making the bosons have 
a common momentum in each polygon. Hence, one as- 
sumes Pi = p,li = I, and <7i = rrii = in Eq.(21), which 
effectively restricts U in Eq.(3) to the pairing- type inter- 
action. Since each K s has its own p and I in Eq.(22), the 
summation must be done. One must sum up all cases by 
replacing each K s in Eq.(22) with J2i P ^s- One therefore 
redefines K s in Eq.(22) as 



l.;p 



( 



\ 



\ 



-u 



i, P 



(23) 

((— l) m in Eq.(18) is absorbed in a negative sign in front 
of [/in Eq.(23).) 

(2) Let us compare the graph like Fig. 4 (Kf 1 ■ ■ ■ Re- 
type) and Fig. 5 (AT„-type). Since each K s of 
K\ x ■ ■ ■ K v s s ■ ■ ■ includes a sum over p and I as in Eq.(23), 
n s (^s) l ' s in Eq.(22) includes the multiple summation. 
Hence, the number of K\ x ■ ■ • A^ s -type graphs is larger 
than that of a K n -type graph. With decreasing tem- 
perature, however, the feature of Bose statistics, accu- 
mulation of particles in the low-energy state, reduces 
these degrees of freedom. At low temperature, all K s 
of K\ x ■ ■ ■ K v s a ■ ■ ■ are likely to have a common p and 
I. Hence, with decreasing temperature, the various 
K? ■ ■ ■ K»° terms (= [£ x(p u h)]* ■••[£ x{PsJs) s ] u ") 
approach a single form K n — J2 P i x (pJ) n i because of 
n = J2 S Sl/ s- 

(3) In an attractive Bose gas (U < 0), K s is always pos- 
itive in Eq.(23), and 2„ in Eq.(22) is therefore a positive- 
term series. Hence, to derive the upper and the lower 
bound for 3 n , one can take only large terms in Eq.(22) 
without considering the cancellation of terms. (For a re- 
pulsive Bose gas, the situation is different [15|.) 
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(a) 



D 2 (2,0) D 2 (0,1) D 3(3,0,0) D 3 (1,1,0) 



a m © 



D 3 (0,0,1) 



FIG. 6: Polygon clusters, (a) For n=2, D 2 (2,0) = 1 and 
Da (0,1) = 1. (b) For rc=3, D 3 (3,0,0) = 1, D 3 (1,1,0) = 1, 
and D 3 (0,0,1) = 3. 



D 4 (4,0,0,0) 



4x 

04(2,1,0,0) 



D 4 (0,2,0,0) 



D 4 (1,0,1,0) 





D 4 (0,0,0,1) 



FIG. 7: Polygon clusters for n=4: D 4 (4, 0,0,0) = 1, 
D 4 (2, 1,0,0) = 2, D 4 (0,2,0,0) = 2, D 4 (l, 0,1,0) = 4, and 
D 4 (0, 0,0,1) = 5. 



B. Graphical consideration of a polygon cluster 

For D n (i/i, . . . , i/ s , . . .) in Eq.(22), it is useful to draw 
explicitly polygon clusters in some cases, and get a feel- 
ing for its magnitude. For n = 2,3,4, we classify the 
distributions {v s } satisfying n — su s , and we obtain 
some D n 's as illustrated in Fig. 6 and 7. 

n = 2. For {^1,^2} satisfying 2 = z/i + 2^ 2 , one has 
{2,0} and {0,1}. 

n = 3. For {^1, ^2, f 3 } satisfying 3 = vi + 2i^ 2 + 3^3, 
one gets {3, 0, 0},{1, 1, 0} and {0, 0, 1}. 

n = 4. For {z^, j/ 2 , ^3, i/ 4 } satisfying 4 = 
t-i + 2// 2 + 3^ + 4^4, one obtains {4, 0, 0, 0}, {2, 1, 0, 0}, 
{0,2,0,0}, {1,0,1,0}, and {0,0,0,1}. 



Fig.6(a) shows D 2 (2, 0) = 1,-D 2 (0, 1) = 1. 

Fig.6(b) shows D 3 (3,0,0) = 1, D 3 (1,1,0) = 1 and 
£> 3 (0,0,1)=3. 

Fig.7 shows D 4 (4, 0,0,0) = 1, D 4 (2, 1,0,0) = 2, 
L> 4 (0,2,0,0) = 2, D 4 (l, 0,1,0) = 4, and D 4 (0, 0,0,1) = 
5. 

In view of Fig. 6 and 7, one notices the following: 
(1) In polygon clusters with n interaction lines, there 
are many types of {v s } ranging from {n, 0, ...,0} to 
{0, ...,0,1}. From now on, we call a group of graphs 
with the same distribution {v s } a "pattern." (For n=4, 
there are five patterns in Fig.7.) The number of patterns 
p(n) appearing in the n th-order graphs is equivalent to 
the number of ways of writing n as a sum of positive in- 
tegers (n = sv s ). In combinatorial mathematics, this 
problem is known as "partition of n" [16| . It is known to 
have the following asymptotic form as n — > 00: 



p(n) 



1 



4V3? 



■ exp 




(24) 




(2) At a given n, the pattern with the largest number 
of configurations D n consists of the largest polygons in 
size (the 2n-gon, {0, . . . , 0, 1}). We can confirm that this 
tendency becomes more evident as n increases. From 
now on, we abbreviate D„(0, . . . , 0, 1) as D(n). (Fig. 4 
and 5 are examples of polygons with seven interaction 
lines. Due to various arrangements of dotted lines, there 
are more Fig.5-type graphs than Fig.4-type ones.) 

Combining p(n) and D(n), one derives an inequality 
for TiD n at a given n in n — ► 00, 

D(n) < ^ D„(vi, ...,v s ,...) 

{".} 



(25) 



where the lower bound is a case in which the sum over 
{vs} in S-D n includes only a pattern composed of the 
largest polygons, and the upper bound is a case in which 
all patterns at a given n are included as if their D n is the 
same as D(n) of the largest polygon. 

(3) With decreasing temperature, large polygons, 
which are successively produced by particle interchanges 
between small polygons, become significant. Because the 
"pairing interaction" becomes significant for Bose statis- 
tics {qi — rrii = in Eq.(21)), directions of arrows in 
such polygons are restricted. As shown in Fig. 5, 6, and 7, 
there appear only special types of graphs, in which every 
two vertices connected by a dotted line are separated by 
odd (not even) number of sides. In combinatorial math- 
ematics, the exact formula of D(n) for such graphs was 
recently found [13], and it has the following asymptotic 
form (Appendix.C) [H[: 



D(n) 



An 



as n 



00. 



(26) 
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Substituting Eq.(26) into Eq.(25) yields an inequality 
for £Z?„ , which will be discussed below. 



C. The instability of an attractive Bose gas in the 
normal phase 

Let us obtain an inequality for Zy(p) in low temper- 
ature. Replacing K" 1 ■ ■ ■ ■ ■ ■ in Eq.(22) by K n (as 
stated in (2) of Sec.3.A), and using Eq.(25) with Eq.(26) 
in Eq.(22), one obtains 




Using Eq.(23) in K n of Eq.(27) and applying it to 
Eq.(20), one obtains [l| 

Z (fx)l[explvf^^-x{p,ir) <Z v (p) 

p,l V n=l U I 

00 exp [■kJ^x) \ 



< ^o(/^) JJexp 
P .i 



V 



^ 16V3n 2 

n— 1 v 



If the infinite series over n in the exponent of the upper 
bound is convergent for all p and I, it guarantees the 
finitcncss of Zy(p), hence the stability of an attractive 
Bose gas. Cauchy-Hadamard's theorem asserts that the 
radius of convergence r c of a power series YlnLi a nX n is 
given by (see Appendix. D) 



lim Kl 1 /" . 



(29) 



Applying this theorem to the upper bound [a n — 

exp(7r^ / /^)/(16v / 3n- 2 )] and the lower bound [a n — 

l/(4n)] in Eq.(28), and using linin^oo n 1 / 71 = 1, one ob- 
tains 1 as r c 's of both bounds. Hence, Zy{p) has the 
radius of convergence r c — 1 as well. For each p, I, the 
convergence condition x(p, I) < r c is given by 



U 



(^-M) 2 +(f)' 



< 1. 



(30) 



The equation of state (1) and (2) implicitly determines 
the chemical potential of a gas as a function of its tem- 
perature and pressure or density. At high temperature, 
/i <C is satisfied for the regions in P — p space, and 
Eq.(30) is therefore satisfied with respect to all p and 
I. Hence, an attractive Bose gas (U < 0) exists as the 
thermodynamically stable state. 

With decreasing temperature, however, the negative p 
gradually increases at a given P or p. Changing P, T 
or p. T in p of Eq.(30) makes possible the violation of 
the convergence condition. Among many p and I, this 



condition is first violated at p 
certain critical value p c as 



I = when \i reaches a 



U 

M 



= 1, 



where p, c satisfies 



= -^-Uk B T c 



(31) 



(32) 



When fi reaches p c , Zy(p) abruptly diverges without 
any precursory behavior. Hence the specific volume v , 
determined by 



1 



knT lim 



d_ 

dp, 



\nZ v (p) \ 
V J 



(33) 



discontinuously jumps to zero. This means that each 
particle cannot keep large distances from other particles. 
Hence an attractive Bose gas undergoes a transition to 
the high-density state. 

The condensation point (P Cl p c ) in the isothermal curve 
is determined by the equation of state and 



p(p c ,T) = -\J-Uk B T 



(34) 



at a given temperature^^. 

For considering the influence of attractive force on the 
coherent many-body wave function, its critical size is of 
great interest. Specifically, the critical size distribution 
h(s) of elementary polygons in Sec.2B is a significant 
quantity [2l[. At high temperature (p ^ — fc^T), h(s) is 
a rapidly decreasing function of s. With decreasing tem- 
perature, h(s) gradually changes to a weakly s-dependent 
distribution, finally reaching h c (s) at the condensation 
line 



Pc 



h c (s) 



X 3 s 



p = 



beyond which distribution an attractive Bose gas does 
not exist in the gas phase. Hence, h c (s) is a critical size- 
distribution of the coherent many-body wave function in 
the gas phase. 

The divergence of Zv(p) occurs first in bosons with 
zero momentum. This means that among many bosons 
with various momentums, the bosons with zero momen- 
tum escape first from a gas, and make a liquid droplet at 
a certain point of the gas. Since this droplet is a high- 
density boson with p = 0, it is likely that the growth of 
this droplet leads to a liquid in the BEC phase. As stated 
in Sec.l, in thermal equilibrium a Bose gas in the BEC 
phase undergoes no GLC, whereas its GLC in the normal 
phase is likely to trigger the formation of the BEC con- 
densate as a by-product. The GLC of a normal helium 
4 gas to a superfluid liquid at P < 0.05 atm falls into 
this category. (One finds a parallel on this point in the 
BCS transition, in which the formation of loosely bound 
pairs of fermions is promptly followed by the formation 
of their BEC condensate.) 
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D. Weak U limit 

In contrast with the classical GLC, the weak-coupling 
limit (U ~ 0) is one of the realistic situations of the 
quantum GLC. Even if —U is very small, the condition 
(34) is satisfied for a small p at the vicinity of the BEC 
transition point in the normal phase. In principle, an ar- 
bitrarily weak attractive interaction creates the quantum 
GLC at sufficiently low temperature. 

When we approximate p in Eq. (34) with po of an ideal 
Bose gas, we can make a rough estimation of T c and p c 
under the fixed density. In an ideal Bose gas, Eq.(2) is 
expanded as [22j 



33/2 ( 



,0^ 



2.612 - 3.545^1/^(7% (35) 



where g a {x) = E„a;™/« a - Using Eq.(34) in p(T) of the 
right-hand side of Eq.(35) at T c , and dividing both sides 
with Xq/v — 2.612 (an ideal Bose gas), one obtains 



s)' +U6 

At U ~ 0, one gets 
and 



k B T \Po 



1 + 0.90 



0.25 



1 = 0. 



-u 

k B T 



Mc ^ -V-Uk B T 



1 + 0.45 



-U 



k B T 



(36) 



(37) 



(38) 



(To denotes the BEC transition temperature of an ideal 
Bose gas at the same p [23|.) 

Figure. 8 shows a schematic phase diagram. The dot- 
ted curve depicts the P — T curve of an ideal Bose gas, 
P = {m/27rh 2 ) 1 - 5 g 5/2 (l)(k B T Q ) 2 - 5 . The solid curve illus- 
trates the P~T C curve (condensation line) derived by ap- 
plying Eq.(37) to the P — T curve. As a crude estimation 
oITbec (traced back to London), we use X 3 /v = 33/2(1) 
with the experimental value of liquid density, and get 
k B T B EC = 2Trh 2 m- 1 2.612- 1 ^{N/V) 1 ^. In Fig.8, the 
solid curve and a straight line T = T BB c intersect at 
the triple point. This schematic diagram is useful only 
for conceptual understanding, and comparing it with the 
phase diagram of the helium 4 needs elaborate calcula- 
tions of some quantities, especially of the liquid state. 



IV. DISCUSSION 

A. Comparison with the classical GLC 

In the phase diagram of helium 4 like Fig.l, with in- 
creasing pressure, the nature of the GLC continuously 
changes from the quantum to the classical GLC [24| . Let 




FIG. 8: A schematic phase diagram. The dotted curve is the 
P — To curve of an ideal Bose gas at the BEC point. The 
solid curve is the condensation line, that is, the P — T c curve 
obtained by applying Eq.(37) to the P — To curve. 



us compare the classical and the quantum GLC. Since 
this paper considers only the attractive interaction for the 
quantum GLC, we can make comparison only on the ini- 
tial stage of the instability. In the classical GLC, Mayer's 
method considers the singularity which must appear in 
the perturbation expansion of Zy(p) = Sa„a; n . In gen- 
eral, the relationship between the coefficient a n and the 
interaction potential is so complicated that it is a difficult 
problem whether the a n derived from a given potential 
actually leads to such a singularity [25| . 

(a) In a classical gas, since the GLC occurs at relatively 
high temperature, the inelastic scattering is important. 
Hence, one must perform an (s — l)-dimensional inte- 
gral for obtaining the cluster integral. In the GLC of a 
dilute Bose gas, however, since it occurs at very low tem- 
perature, the diagram with a common momentum p like 
Fig. 5 becomes important. Hence, the integrals of a 2s- 
sized polygon are factorized, and the bounds for Zy{p) 
become simple as in Eq.(28). 

(b) In a classical gas, one can prove the convergence 
or the divergence of Zy(p) only after summing up all 
momentum components. In an attractive Bose gas, how- 
ever, the instability is expected to occur at bosons with 
zero momentum. Hence, one can prove the divergence of 
Zy (p) by focusing on the p = component of 3„ . 

(c) In the classical GLC, due to the absence of Bose 
statistics, there are far more Kf 1 ■ ■ ■ K^'-type clusters 
than K n -type polygons. Hence, at each order of the per- 
turbation, one cannot estimate Zy(p) only by the largest 
polygons in size. For the quantum GLC, however, one 
can estimate Zy(p) by focusing the largest polygons as 
in Sec.3.B. 

These features (a) ~ (c) will become useful for eluci- 
dating the essence of the quantum GLC and studying the 
classical-quantum crossover of the GLC in future studies. 
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B. The quantum GLC in fermions 



There is another type of quantum GLC occurring in 
fermions such as a helium 3 gas. In the phase diagram of 
helium 3, the location of the normal-superfLuid transition 
in the liquid phase is far apart from the condensation line 
between the liquid and the gas phase. (This is in contrast 
with helium 4 in Fig.l.) Hence, it is uncertain to what 
extent the observed GLC is affected by Fermi statistics 
0. In principle, however, the quantum GLC in fermions 
has some different features from that in bosons. For an 
ideal Fermi gas, due to one-particle excitation, the anti- 
symmetrical coherent wave functions are not formed even 
at zero temperature, whereas it shows the Cooper insta- 
bility with the attractive force. Hence, it is appropriate 
to begin with the pairing interaction when considering 
the GLC. In a Fermi similar argument to that in 

Sec.2.C is valid as well: When two Fermi particles (p 
and p') belonging to two different bubbles have the same 
momentum (p — p'), a new diagram containing a square 
must be included with a negative sign [26| [27]]. Its K s 
has a similar form to Eq.(23) except for an odd integer I, 
and 5 n of an attractive Fermi gas has a similar structure 
to Eq.(22) in appearance. But the argument in Sec. 3. A 
is affected by Fermi statistics as follows. 

(a) Due to Fermi statistics, each K s in Eq.(22) has 
its own momentum and frequency even at zero temper- 
ature, and the number of the K" 1 ■ ■ ■ K^-type cluster is 
therefore larger than that of the K n -type polygon in H„. 
Hence, at a given n, dominant terms in the perturbation 
series do not come from the K n -type polygon, but from 
the various Kf 1 ■ ■ ■ K^-type polygon clusters. (In n=7 
for example, the polygon cluster like Fig. 4 is more impor- 
tant than the single largest polygon like Fig. 5.) In this 
respect, the GLC in a Fermi gas has the same feature as 
the classical GLC [see (c) in Sec. 4. A]. 

(b) For Fermi statistics, a negative sign appears in 
front of U in K s (Eq.(21)), because of a) p a*_ p a^ p ia p > = 

— aJa_p'aL p Op' for p ^ p' . Accordingly, the negative sign 
in front of U (< 0) in Eq.(23) disappears. Hence, H n of 
an attractive Fermi gas is an alternating series. The 
perturbation series of an attractive Fermi gas has a com- 
pletely different structure from that of an attractive Bose 
gas. That is, one must consider the cancellation of terms 
which are opposite in sign, and cannot approximate 5 n 
only by the terms with large D n . One must carefully 
count all K" 1 ■ ■ ■ K" a -type clusters in Eq.(22). Gaudin 
[2(| and Langer [271 ] gave an approximation of Zy([i) 
along this line of thought, and led the BCS gap equation 
from it. (Using their result, a rigorous proof was given re- 
cently that, as long as attractive force is increased within 
the BCS scheme, the GLC is impossible 28].) Applica- 
tion of the inequality method to an attractive Fermi gas 
is a future problem. 



C. Future problems 

Lastly, we point out some problems in future studies 
of the quantum GL C . 

(a) For Mayer's method of condensation, we know a 
long-standing problem in the interpretation of results. 
Since this method focuses on the instability of a gas and 
contains no information on a liquid, there remains a ques- 
tion whether the singularity of the power series in Zy (/i) 
actually gives the condensation temperature T c , or corre- 
sponds to an end point of supersaturation T s . This prob- 
lem is related with technical difficulties in estimating the 
grand partition function. In general, it is difficult to treat 
the highly inhomogeneous system so that it will reach 
true thermal equilibrium. Poor approximations of Zy(fi) 
often include the uniform-density restraint, thus prevent- 
ing the system from reaching true thermal-equilibrium. 
The singularity of Zy{jj) obtained by such an approxi- 
mation corresponds to the end point of supersaturation. 
If a satisfactory approximation is made, it will lead to 
the condensation point in thermal equilibrium. Hence, 
the validity of Mayer's method depends on the type of 
phenomenon. For bosons, we can have the following op- 
timistic view: The uniformity of the system by p = 
particles is the characteristic property of bosons at low 
temperature. Hence the satisfactory approximation of 
Zy(fi) is more probable than that of more inhomoge- 
neous systems. A quantitative analysis of this point is a 
future problem. In this paper, we used the term "T c of 
the GLC" as a generic name for the temperature at which 
a Bose gas shows the instability leading to the collapse. 

(b) To discuss a liquid on a common basis with a gas, 
one must make a more realistic model by including the 
short-range repulsive interaction into Eq. (3) . In this case, 
one must expect interference of the attractive U and the 
repulsive U' in 3 n , with a result that 3 n becomes inter- 
mediate between the positive-term series and the alter- 
nating series. An estimation of such a series is a difficult 
problem. 

(c) For the classical GLC, there appears long-range cor- 
relation at the vicinity of the condensation point, which 
comes from the general mechanism of a first-order phase 
transition. For the quantum GLC in an attractive Bose 
gas, the long-range coherence comes from Bose statistics. 
This difference may affect the nature of fluctuations in 
the vicinity of a transition point. 

(d) The recent experiment of trapped atomic Bose gas 
will provide us with a new tool (29|. Ultralow tempera- 
ture is available and direct control of the effective inter- 
particle force becomes possible using magnetic field. The 
system of interest from our viewpoint is not a BEC gas, 
but a normal gas just above the BEC transition point 
[30| . The trapped atomic Bose gas, however, has its 
own properties which may complicate the foregoing argu- 
ments, that is, metastability. The confinement of atoms 
by a harmonic potential gives rise to the spatial inhomo- 
geneity in the system. Many bosons gather at the cen- 
ter of a trap. The resulting increased density enhances 
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the zero-point motion of particles, which stabilizes the 
gas state. For a trapped atomic Bose gas with attrac- 
tive interaction, one finds a metastable BEC state at low 
temperature and below critical density. Up to now, the 
collapse of such a metastable BEC gas was detected [3lJ 
but the GLC in the normal phase has not been ob- 
served. Further, there is a possibility that a state after 
the transition is not a liquid but a solid or molecules. If 
it is true, we must consider this system from a different 
viewpoint. 
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APPENDIX A: COMMENTS ON REF.[8] 

As a starting point to Zy(u) of an attractive Bose gas, 
Ref. [8j began with an approximate formula [Eq.(9) in 
Ref. [gj], which carefully deals with the Kf 1 ■ ■ ■ K%"-type 
cluster like Fig4. As confirmed in Sec. 3, however, the pe- 
culiar nature of a Bose gas appears in the n — > oo limit. 
In this limit, Eq.(9) and all subsequent integrals (with 
respect to t) in Ref.[j| turn out to diverge. The finial 
form of Zy(u) (Eq.(21) in Ref.Q) approaches oo — oo, 
which loses the mathematical meaning. Hence, it is ques- 
tionable to identify it as a Yang-Lee zero. (One way to 
avoid this divergence is to perform all integrals in Ref. 8] 
over [0, M] with a finite M instead of [0, oo]. But such a 
trick is unnatural in a Bose gas.) In this paper, we prove 
the singularity of Zy (//) by showing its divergence using 
the inequality; this gives a sound basis to the instability 
condition (Eq.(34)), which was first obtained in Ref.0]. 



APPENDIX B: POLYGONS IN Z Q (jt) 



After extracting T(p = 0) from the integrand of Eq.(B2) 
and performing the integral, one obtains 

f.= (x* + A.\X>- l >Z s ). (B3) 



APPENDIX C: ASYMPTOTIC FORMULA OF D 

The asymptotic form of D(ri) (Eq.(26)) is derived as 
follows [L7T ] . We think of regular 2n-polygons with n in- 
teraction lines, in which n lines connect n pairs of vertices 
separated by an odd number of sides. We denote its num- 
ber by D(n), and classify it by a transformation property 
such as rotation and reflection: A polygon is invariant 
under rotations through angle ixljn (I = 0, 1, . . . , In — 1). 
Further, a 2n-polygon is symmetrical with respect to the 
2n axis. Hence the symmetry group is composed of An el- 
ements Co, C\, . . . , C*4 n _i. We denote by 0(Cj) the num- 
ber of graphs which are invariant under transformation 
Ci. 

Birnside's theorem [H[ states D(n) = ^(d)/4n. The 
most important </>(Cj) is the number of polygons invariant 
under an identical transformation 4>(Cq) = n\ (the num- 
ber of ways of connecting n pairs of vertices). Among 
many </>(Cj)'s, the 4>(Cq) becomes dominant as n — > oo. 
Hence, taking only 4>(Cq) in E0(C,-), one obtains an 
asymptotic form D(n) — > n\/An at n — > oo. 

For D(n) with a finite n, one must classify 2rt-polygons 
with interaction lines by symmetry, and include other 
0(Ci)'s in E<f>(Ci). The result is as follows [13]. 

Theorem 1 (Nakamura) The number of the nonequiv- 
alent 1-regular odd spanning subgraphs of the complete 
graph of order 2n by the Dihedral group £>2n is given by 

I i 2n ~ 1 i 
D ^ = Tn + -n^^ + -^ + S ^ 

4=1 

where 



f s in Eq. (9) is estimated as follows [Il| . If x 2 sl [the last 
term in the exponent of Eq.(9)] is replaced by x 2 , f s be- 
comes a function of x\ — Xo(= xiq) as / s (^io)- Equation 
(9) shows that f s (xio) is made by (s — 1) times of the 
convolution of exp(—mx 2 /2(3h 2 ). Hence, using its three- 
dimensional Fourier transform T(p) = A' 3 exp(— X' 2 p 2 /2), 
where A' = X/\^2tt, f s (x\o) is expressed as 



fs(x 10 ) 



V 



2tt) ; 



■ J (V^^^TipYe-^^A^dp, 



(Bl) 

where V corresponds to the translation of the elementary 
polygons. Since f s — f s (xio = 0), one has 



f s = V (V2ny s T(p) 



. 4irp 2 dp 
(2tt) 3 ' 



(B2) 



r o, 



R 2n = { 



^ 2 s sW. 

d=2s+t 



G)" 



d = odd, n = 2m 



d = odd, n = 2m + 1 



2n 



d/2 



5 '■ 



d = even, 



(d is the greatest common divisor of 2n and i, and s 
and t are positive). 



0, 



n = 2m 



2 (n-i)/2(n^ ]{ n = 2m +l, 
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2 s sW. 



n=2s+t 



APPENDIX D: CAUCHY-HADAMARD'S 
THEOREM 

We begin with a general series J2^Li c n- If one can 
find k and m such that ^/cvT < k < 1 when n > m, one 
gets for a large n 



Since < k < 1, S^Li c " converges. 

The above result is used in XmLi flni n with |a n a;™| 
identified as c n . By defining linin^oo la,^ 1 /™ = I, one 
obtains linin^oo lanx"! 1 /™ = l\x\, that is, lim^-joo ^/c^" = 

z|z|. 

This means that, if l\x\ < 1, Eq.(Dl) is valid, hence 
Y^=i°n converges. Contrarily, if l\x\ > 1, it diverges. 
The radius of convergence r c of X^^Li a nX n is given by 
l\r r \ — 1, that is. 



hm \a n \ 1 / n . 



(D2) 



>S n = Ci H h C„ < fc + 



+ fc"= fc( ; P- . (Dl) 
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